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Abstract 

We present a detailed study of prompt photon production cross section in heavy- 
ion collisions in the central rapidity region at energy of y/s = 5.5 TeV, appropri- 
ate to LHC experiment. We include the next-to-leading order radiative corrections, 
0{a em a s ), nuclear shadowing and the parton energy loss effects. We find that the 
nuclear effects can reduce the invariant cross section for prompt photon production 
by an order of magnitude at pt = 3 GeV. We discuss theoretical uncertainties due to 
parton energy loss and nuclear shadowing parameters. We show that the K-factor, 
which signifies the importance of next-to-leading order corrections, is large and has 
a strong p t dependence. 



1 Introduction 



There has been a considerable theoretical and experimental interest in studying photon 
production in heavy-ion collisions at BNL's Relativistic Heavy Ion Collider (RHIC) and 
CERN's Large Hadron Collider (LHC) energies [I]. Photons produced in heavy-ion colli- 
sions provide an excellent probe of the properties of the dense matter, such as the quark- 
gluon plasma or the hot hadronic gas, produced after the collision. Due to the small cross 
sections of electromagnetic interactions, photons can escape the strongly interacting matter 
produced in the collision without further interactions. 

Studying photon production at the RHIC and LHC energies is of special interest, as 
it has been suggested as an elegant signal for detecting the formation of a quark gluon 
plasma (QGP) in heavy-ion collisions 0. However, photons can be produced at different 
stages of the heavy ion collision and thus have different origin. For example, photons can 
be produced at the early stages of the collision through QCD processes such as qg — > qj 
or they can be emitted from a thermalized quark gluon plasma or hadronic gas. Another 
source of produced photons is decay of hadrons such as pions and etas produced in the 
heavy ion collision ||. Furthermore, different processes give the dominant contribution at 
different p t 's. Therefore, it is quite difficult to make reliable predictions for the absolute 
number of photons produced in a heavy ion collision [fj. H. 

Prompt photons are an important background to thermal photons in the low to in- 
termediate pt region and are dominant in the high p t region. Therefore, it is extremely 
important to be able to calculate their production cross section reliably. Fortunately, one 
can use perturbative QCD in the high p t region to calculate the prompt photon production 
cross section. In this work we continue our study of the production of high p t (p t > 3 GeV) 
prompt photons in heavy ion collisions Q. Prompt photons are produced either directly in 
the hard collisions of the partons inside the nuclei like qg — > q^ or through bremsstrahlung 
of quarks and gluons produced in the hard collision such as qg — > qgj. We include all next- 
to-leading order, 0{a em a 2 s ) QCD processes as well as nuclear shadowing and medium 
induced parton energy loss effects. We discuss and estimate the theoretical uncertain- 
ties due to different choices of nuclear shadowing and energy loss parameters. We show 
that next-to-leading (NLO) corrections are large and must be included to make reliable 
predictions. 

In section I we review the prompt photon production in hadronic collisions in next- 
to-leading order. Based on QCD factorization theorems, we write down the expression 
for prompt photon production and list the hard partonic processes involved. In section II 
we discuss the nuclear effects such as shadowing and energy loss involved in production of 
prompt photons in heavy ion collisions. In section III we present our results for the prompt 
photon production invariant cross section, E-jj§-, at LHC energies and an estimate of the 
theoretical uncertainties due to variation of nuclear shadowing and energy loss parameters. 
We show that the effective K-factor, defined as the ratio of NLO to LO cross sections 
in heavy ion collisions is large and has a strong p t dependence. This clearly shows the 
importance of including next-to-leading order contributions to prompt photon production. 
We conclude with a discussion of prospects for detecting nuclear effects by measuring 
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prompt photons at LHC. 



2 Prompt photon production in pQCD 

Using factorization theorems and perturbative QCD, the inclusive cross section for prompt 
photon production in a hadronic collision can be written as a convolution of parton densities 
in a hadron with the hard scattering cross section and the parton to photon fragmentation 
function: 



£ 7 J^(v^P 7 ) = jdx a jdx b jdz ]T Fi(x a , Q 2 )Fj(x b , Q 2 )D ch (z, Q))E, °^ cX (1) 

where the Fi(x,Q 2 ) is the i-th parton distribution in a nucleon, x a and x b are the frac- 
tional momenta of incoming partons, D c / 1 {z ) Q 2 ) is the photon fragmentation function 
with z being the fraction of parton energy carried by the photon. The parton-parton cross 
sections, d G ^ cX , include all processes up to and including 0(a em a^), such as leading-order 
subpro cesses: 

q + q -> 7 + 3 

q + g -> i + q (2) 

and the next-to-leading order subprocesses, 

q + q — > g + g + 7 

q + q — > g + g + 7 

g + g' — > g + g' + 7 

g + g — > g' + g' + 7 

g + g' — > g + g' + 7 (3) 

We refer the reader to J?J for a complete list of all 0(a em a 2 s ) processes. There are two 
classes of subprocesses, "direct" production, which does not have convolution with the frag- 
mentation function and the "bremsstrahlung" contribution that has convolution with the 
fragmentation function. Direct subprocesses contribute to leading order as well as next-to- 
leading order, while bremsstrahlung processes only contribute at the next-to-leading order. 
Clearly, only the sum of these two contributions constitutes the full NLO calculation of 
the direct photon production. It is important to note that in order to include all 0(a em a 2 ) 
processes, it is necessary to include processes which formally look order 0{a em a\). This 
is because some of the terms in (§) have a divergence proportional to l/a s which makes 
those processes 0(a em a s ). This divergence comes from the quark and emitted photon 
being collinear. The reader is referred to || for a discussion of the collinear divergences. 

The nucleon structure functions, Fi(x, Q 2 ), and parton to photon fragmentation func- 
tions, D c / 1 (z,Q'j), are also evaluated at the next-to-leading order. We use the MRS99 
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set for nucleoli structure functions and Bourhis et al. parameterization of the photon 
fragmentation functions ||10|| . Structure functions, fragmentation functions and the run- 



ning coupling constant depend on the factorization, fragmentation and renormalization 
scales respectively which are usually taken to be the same and proportional to the photon 



transverse momentum p t . Aurenche et al. [11 have studied the dependence of the prompt 



photon cross section on the choice of scale and have shown that the choice of Q = p t /2 gives 
a very good description of prompt photon production in hadronic collisions. Therefore, we 
will use Q = p t /2 in our calculation. The running coupling constant a s (Q 2 ), calculated to 
next-to-leading order, is given by 



2 12tt r 6(153- 19A^ / )lnlng 2 /A 2 

®s\l% ) ~rrrz - ~ . . ~ 1 r~rr « - ~ . . « 



(33 - 2N f ) In Q 2 /A 2 L (33 - 2N f ) 2 In Q 2 /A 
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(4) 



where Q 2 is the renormalization scale, A is the QCD scale parameter and Nf is the number 
of flavors. 

In Figure (|T]a) we illustrate the importance of including the next-to-leading order con- 
tributions by calculating the ratio of NLO to LO cross sections, the so-called K-factor. 
Clearly, the NLO corrections are large and p t dependent. It would be useful to go even 
beyond the NLO to make sure that the higher order corrections are not even larger. It is, 
however, clear that a leading order calculation in the LHC kinematic region is meaningless 
since the next-to-leading order corrections are huge. 

In Figure (|l|b) we show an alternative definition of the K-factor defined as K = 
NLO/(LO + Brems.) sometimes used in the literature M. The main reason for this def- 
inition was inclusion of an incomplete set of Bremsstrahlung diagrams in previous works 
on prompt photon cross sections |4J] and we show it here for comparison. 

In principle, in the region where x% ~ Pt/ \fs is large, i.e. in the very low and very high 
p t region of phase space, one would need to resum certain logarithmic terms, of the form 
a s ln(l — Xt) and a s \n.x t - However, these resummations are not very important for the 



range of p t and ^fs that are considered in this work and will be neglected |12 

There is also an alternative approach to calculating prompt photon cross sections which 
uses the LO cross sections as well as LO structure and fragmentation functions. In order to 
reproduce the experimental data in this approach, it is necessary to include a phenomeno- 
logical parameter which is loosely identified as the "intrinsic" momentum of the initial state 
partons |13|]. This intrinsic momentum is generated by the initial state radiation of quarks 
and gluons. There is however no theoretical calculation of these effects at the moment and 
one has to model them by generalizing the standard definitions of the parton distribution 
functions from q(x, Q 2 ) and G(x, Q 2 ) to q(x, k 2 , Q 2 ) and G(x, k 2 , Q 2 ). One then introduces 
a weight function, typically a Gaussian, with a width < k 2 > which represents the intrinsic 
transverse momentum of the initial state partons. Assuming some reasonable range in kt, 
these intrinsic momenta are integrated over in the final result. This approach gives fairly 
good description of most of the experimental results. 

It is important to realize that this intrinsic momentum kt grows with energy and can 
be as large as ~ 1 — 2 GeV in fixed target experiments and ~ 5 GeV in the Tevatron WA . 



Therefore, strictly speaking, it is not an intrinsic momentum. In addition, the value of k t is 
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Figure 1: a) The hadronic K-factor defined as the ratio of NLO prompt photon cross 
section, E-^ to the LO cross section, and b) the K-factor defined as the ratio of NLO to 
the LO plus bremsstrahlung cross sections. 



process dependent and thus can not be taken as universal parameter. It might be possible to 
extract its value from the data in different experiments, but its value can not be predicted. 
In our work, we will not follow this approach but instead we will use NLO perturbative QCD 
formalism because it seems theoretically more rigorous and self-consistent. We note that 
NLO calculations give fairly good agreement with the experimental data for p t > 3.5GeV 
at lower energies (y/s < 63 GeV) [[□]] and for p t > lOGeV at higer energies (y/s = 1800 
GeV) ||14|| . However, it is worth pointing out that there are claims of inconsistencies among 
various data sets from different experiments at low energies [O, TA\ . In addition, one should 
keep in mind that our study extends down to lower values of pt, where NLO predictions 
may not be reliable. 



3 Prompt photon production in heavy ion collisions 

To calculate the prompt photon production cross section in heavy ion collisions, we will 
use eq. ([I]) modified for nuclear effects: 



E.^—i^pt) = j dx a J dx b J dz g F t A (x a , Q 2 )Ff(x b , Q 2 )D ch (z, Q*)E y ^ cX (5) 

where the F^(x, Q 2 ) is the i-th parton distribution in a nucleus and D c / 7 (z,Q 2 ) is the 
photon fragmentation function in a nuclear environment. The partonic processes (a partial 
list) are given by (M |3|) as before. The nuclear structure function F^(x, Q 2 ) and fragmenta- 
tion function D c /^(z, Q 2 ) are not known to next-to-leading order, but rather they are NLO 
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nucleoli structure and fragmentation functions modified for nuclear effects. Strictly speak- 
ing, therefore, our calculation of prompt photon production in heavy ion collisions is not 
a next-to-leading order calculation but is the most complete calculation performed so far. 
Below, we will describe the nuclear modifications to structure functions and fragmentation 
functions in detail. 



3.1 The Nuclear Shadowing Effect 

Calculation of the prompt photons production cross section in nuclear collisions requires 
knowledge of the nuclear structure functions F^(x, Q 2 ), usually measured in Deep Inelas- 
tic Scattering (DIS) of electrons on nuclei. It is an experimental fact that Ff-(x, Q 2 ) ^ 
AF { N (x, Q 2 ) where F t N (x, Q 2 ) is the free nucleon structure function. This modification has 
a strong x dependence and is due to having different nuclear effects in different region of 
phase space. At small values of x, for instance x <~ 0.07, the nuclear structure function 
is less than nucleon structure function scaled by A. This is known as shadowing. As x 
grows bigger, nuclear structure functions get bigger than the free nucleon structure func- 
tion. This is known as anti-shadowing. As x further increases, nuclear structure functions 
become less than the free nucleon ones again which is known as the EMC effect. In this 
section, we will concentrate on shadowing and anti-shadowing regions since those are the 
kinematic regions where the prompt photon production takes place. We refer the reader 



to [[L5[] for a recent review of nuclear shadowing. 

An intuitive explanation of nuclear shadowing and anti-shadowing effects in DIS de- 
pends on the reference frame of the nucleus. Of course, the physical observables such as 
structure functions can not depend on our choice of reference frame. However, working in 
different frames helps one identify the different physical mechanisms involved. In the rest 
frame of the nucleus and in perturbative QCD, shadowing of nuclear structure functions 
in DIS is due to multiple interaction of the qq component of the photon wave function 
(emitted by the electron) with the nucleus. The amplitude of qqA interaction is mostly 
imaginary at small x and multiple interactions of the pair with the nucleus introduces a 
phase difference between different amplitudes which leads to a destructive interference. 
This in turn reduces the nuclear cross section. In a non-perturbative description of shad- 
owing, the photon is resolved in terms of its hadronic fluctuations which in turn multiply 
interact with the nucleus. The multiple interactions again reduce the nuclear cross sections 
due to destructive interference. In this frame, anti-shadowing at larger x is due to a large 
real part of the interaction amplitudes which interfere constructively with the imaginary 
part and lead to an enhancement of the nuclear cross section (structure functions). 

In the infinite momentum frame where the nucleus is moving very fast, shadowing is 
caused by high parton density effects small x. The small x partons have a large longitudinal 
wavelength and can spatially overlap and recombine. These recombination effects reduce 
the nuclear parton number densities and hence the nuclear cross sections. Working in this 
frame enables one to treat nuclear shadowing and parton saturation in nucleons on the 
same footing due to the identical physical mechanism involved in both. Anti-shadowing is 
due to longitudinal momentum conservation (momentum sum rule) in this frame. 
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Even though there has been considerable amount of theoretical work done on nuclear 
shadowing and impressive progress made in understanding the physical principles of nuclear 
shadowing 
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we are far from having a precise and quantitative description of nuclear 
shadowing. In practice, one measures the nuclear structure functions in deep inelastic scat- 
tering of leptons on nuclei [HJ. The measured structure functions are then used in nuclear 
collisions. The scale dependence of the nuclear structure functions is even less understood 
due to the limited range of Q 2 covered in fixed target experiments. Also, shadowing of 
gluons is not well understood due to the fact that they can not be directly measured in 
DIS experiments. The working assumption is that high parton density effects are negligible 
and DGLAP evolution equations are valid in which case the gluon distribution function 
can obtained from the scaling violation of the F2 structure functions. This assumption, 
however, will break down at small values of x due to high parton density effects |l7j and 
one will need to measure the gluon distribution function differently. 

In this work we will use two different parameterizations of nuclear structure functions 
due to Benesh, Qiu, Vary [|T^] and Eskola et al. |TJ|]. Both parameterizations fit the current 
experimental data quite well even though they are somewhat different. Also, since there are 
no experimental data at the small x, high Q 2 (Q 2 > lGeV 2 ) region, any parameterization of 
nuclear structure functions in this region is subject to large uncertainties. This is somewhat 
important for RHIC but becomes crucial for LHC. An eA collider such as the proposed 
eRHIC is urgently needed and would greatly improve our knowledge of nuclear structure 
functions in the small x, high Q 2 region as well as reducing the theoretical uncertainties in 
the larger x region. 

The parametrization of the nuclear shadowing function proposed by Benesh, Qiu and 
Vary is given by 



1* 



S(x,A) 




014X 

n ™ N 1+V* 2 (1/ 

- a^oJi-r 



-i/go) 



+k q A a i(l/x-l/x Q ) 



x < x < 0.6 

X < Xq 



(6) 



It gives a good description of all EMC, NMC and E665 data [16]. The parameters k, 



"q, 



aii, ct 2 , a 3 and x are fitted to deep inelastic data for the ratio F 2 A (x,Q 2 )/F 2 D (x,Q 2 ) and 



can be found in p8[ . In this parametrization, nuclear structure functions are independent 
of Q 2 and shadowing of gluons is assumed to be the same as quarks. A more recent 
parameterization of the nuclear structure function is that of Eskola et al. |19| that also fits 
the existing experimental data quite well, has Q 2 dependence and distinguishes between 
quark and gluon structure functions |Tj|. We show the nuclear shadowing ratio defined 
as S = F^/AF^ in Figure (Q) using the parameterizations of BQV [|18| and EKS ]19 . 
Clearly, the two are quite different although both fit the experimental data fairly well. 
This signifies the need for a high energy collider experiment such as eRHIC to improve on 
the current measurements of the nuclear structure functions. 

In this work, we will use both parameterizations and investigate the dependence of 
prompt photon production cross sections on our choice of nuclear structure functions. 
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Figure 2: The nuclear shadowing ratio as parameterized by BQV [|L8] and EKS98 [|19|] ). 
3.2 Medium induced parton energy loss effects 

Another difference between hadronic and heavy ion collisions is in the multiplicity of the 
final state particles produced. In high energy heavy ion collisions and in central rapidities, 
the multiplicity of particles produced per unit rapidity is much larger than that in hadronic 
collisions. Therefore, many body effects such as secondary collisions, which are negligible in 
hadronic collisions become important in high energy heavy ion collisions. Another example 
is the medium induced energy loss. In hadronic collisions, since the number density of 
particles per rapidity produced is small, one can neglect the multiple interactions of the 
produced partons with each other. On the other hand, in high energy heavy ion collisions, 
due to high multiplicities per unit rapidity, the multiple interactions of produced partons 
with each other can not be neglected. As a consequence of multiple interactions with 
the medium, produced particles can lose energy before hadronizing. This affects their 
energy and momentum spectrum. Energy loss of energetic partons passing through a 
dense medium has been a hot topic lately |Z0 |. There has been a considerable progress 



made in understanding the different forms of the energy loss in different limits. It has been 
shown, for a finite size medium, that parton energy loss increases with increasing parton 



energy f2"0|| . Current calculations of energy loss effects are done at the leading order 0(a s ) 
and a next-to-leading order calculation is not presently available. 

A rigorous treatment of energy loss effects in a heavy ion collision is extremely com- 
plicated and beyond the scope of this work. Rather, we will take a phenomenological 
approach to medium induced energy loss in nuclear collisions and use a model of Wang, 



Huang and Sarcevic [pi], to estimate energy loss effects in high energy heavy ion collisions. 
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In this commonly used model, it is assumed that the main effect of multiple interactions 
in the medium can be accommodated by modifying the photon fragmentation functions. 
In the central rapidity region, parton produced in the hard collision is traversing the nu- 
clear medium and losing its energy as a result of multiple interactions with the deconfined 
medium. This parton will hadronize outside the nuclear medium but with a reduced energy. 

In the energy loss model of Wang, Huang and Sarcevic pi| , the parton to photon 
fragmentation function, zD®,(z, Q 2 f), which gives the probability for a parton to fragment 
into a photon, is modified to include multiple scatterings of the fragmenting parton from 



the nuclear medium before it fragments. The nuclear fragmentation function zD a / 1 (z } Q 2 
is given in terms of the photon fragmentation function zD®,(z, Q 2 ) by [21 



zD a/7 (z,AL,Q 2 f ) 



1 N 

N n=0 



n , 



t a n Dl h (z:,Q}) + J24D g/l (4,Q}) 



3=1 



(7) 



where z® = z/(l — (J2i=o e i) / Pt) ■> z$ = zE t /e°j and P a (n) is the probability that a parton 
of flavor a traveling a distance AL in the nuclear medium will scatter n times. It is given 
by 



Pain) = !^e-^, (8) 
n! 

The first term in Eq. (|7|) corresponds to the fragmentation of the leading parton a with 
reduced energy E t — Yh=o e a after n gluon emissions and the second term comes from the 
j-th emitted gluon having energy e J a , where is the energy loss of the parton a after j-th 
scattering. Since we are studying the energy loss effects only phenomenologically, we will 
consider two different cases for the energy loss per collision, e J a . First we will take it to 



be constant, as considered for example, in |22]. We will also consider energy dependent 



energy loss in the form of e{ = a s yfi 2 \ a El, where E{ is the energy of the parton a after 
j scatterings, A a is the inelastic mean free path of parton a and fi 2 represent a screening 
mass generated by the plasma and serves as an infrared cut off. It should be emphasized 
that the general form of the energy loss per scattering, e J a , is theoretically unknown and we 
consider two possible cases. In order to study theoretical uncertainties involved, we will 
take several values for the constant energy loss and in the case of energy-dependent energy 
loss we will vary parameters /i 2 and A a . 

The nuclear fragmentation functions of quarks and gluons as defined in (|7|), with en- 
ergy dependent energy loss, are shown in Figure Here we have taken /x 2 = lGeV 2 
and A a = lfm. The constant energy loss case is quite similar. For the photon fragmen- 



tation function, zD®,(z, Q 2 ), we use the parameterization of JlO|. The parton to photon 
fragmentation functions in a nuclear medium are enhanced at small z and suppressed at 
large z as compared to the fragmentation functions in vacuum. This is due to the fact 
that high energy (high z) partons multiply scatter from the medium and lose their energy 
which shifts their energy fraction z to a smaller value. The nuclear fragmentation function 
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obtained assuming constant energy loss per scattering has qualitatively the same behavior 
function of z. 




Figure 3: The photon fragmentation function, a) zD®,^(z, Q 2 ) (solid line) fT(| and 
the nuclear fragmentation function zD q / 1 (z,Q'j) obtained using (|7|) (dashed line) and 
b) zD®,(z, Q 2 ) (solid line) |L0| and the corresponding nuclear fragmentation function 
zD g /^(z,Q 2 ) (dashed line). 

As one goes to higher energies in heavy ion collisions, one probes smaller and smaller en- 
ergy fractions z in the fragmentation functions. The fragmentation functions, zD®^(z, Q 2 f), 
are fit to the experimental data and parameterized. However, the existing data does not 
cover the z ranges which will be explored by LHC. Therefore, the current parameteriza- 
tions of parton to photon fragmentation functions are set equal to zero below some energy 
fraction, z < Zq ~ 0.01. This is shown in Figure (|3|). In the kinematic region appropriate 
to LHC, however, on will need to know the fragmentation functions below this energy 
fraction. Therefore, we use both the standard (hadronic) fragmentation functions which 
are zero below zq and another parameterization which is identical to the standard one for 
z > z but is set equal to a constant for z < zq, i.e. zD^ g ^(z,Q 2 ) = zD® g ^(z = zq,Q 2 ) 
for z < Zq. We have checked that the difference between the two parameterizations leads to 
a negligible (< 1%) difference in the nuclear cross section. The reason is that the average 
z is sufficiently large (z ~ 0.1) and thus the results are not sensitive to the variation of 
fragmentation functions in the small z region. Also, partonic cross sections at small z are 
power suppressed and do not contribute significantly. 

The parton mean free path A a and the screening mass fi 2 are largely unknown and 
in principle depend on the parton species and the medium temperature. In this work, 
we treat /x 2 and A a as unknown parameters and show the dependence of our results on a 
physically reasonable variation of them in Figure (f|). By varying these parameters, fi 2 and 
X a , we study theoretical uncertainties. Experimental determination of these parameters 

would be difficult. Energy dependent form of energy loss, e{ = a s \J {i 2 \ a E[, is, strictly 
speaking, valid only for coherent photon radiation. Here, we have considered this form just 
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Figure 4: The uncertainty in prompt photon cross section at ^/s = 5.5 TeV due to variation 
of energy-dependent energy loss parameters. 

as a phenomenological expression that gives photon p t distribution that we could compare 
with the constant energy loss case. The precise form of energy loss in a realistic nuclear 
collision, such as those at RHIC and LHC energies, is not well understood at the moment 
and is expected to be extremely complicated. 



4 Discussion 

We show our results for the prompt photon cross section in nuclear collisions at LHC 
energies in Figures @, @ and ([/]). We use different forms of nuclear shadowing and 
energy-dependent energy loss with fi 2 = lGeV 2 and A a = lfm (Figures 5 and 6) or a 
constant energy loss (Figure 7). We find nuclear effects at LHC to be strikingly large. 
The nuclear cross sections can be reduced by 90% at p t = 3 GeV and 50% at p t = 20 
GeV for the energy dependent energy loss scenario. The constant energy loss of 2 GeV per 
scattering leads to the similar size effects while smaller energy losses per scattering, such 
as 1 GeV, 0.5 GeV and 0.25 GeV lead to smaller suppression of the nuclear cross section. 
In Figure (^) BQV parameterization of shadowing is used which is Q 2 independent. This 
reflects itself in almost constant (~ 30%) contribution of shadowing to the suppression of 
the nuclear cross section at all pt- On the other hand, energy loss effects become smaller 
at larger p t 's (~ 50% at p t = 3 GeV and ~ 20% at p t = 20 GeV). 

In Figure (|]) we use both the BQV JTjJ and the EKS98 [HJ parameterization of shad- 
owing. The EKS98 parameterization has Q 2 dependence and distinguishes between quark 
and gluon shadowing while the BQV parameterization is Q 2 independent and does not 
distinguish between quark and gluon shadowing. The difference at low p t (pt ~ 3 GeV) 
is small, about 4% while at p t = 20 GeV it is about 20%. This is mainly due to the Q 2 
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dependence of EKS98 parameterization. 
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Figure 5: Prompt photon cross section in the central rapidity region at \fs = 5.5 TeV 
without nuclear effects (solid line), with BQV shadowing and no parton energy loss 



(dotted line), with parton energy loss and no shadowing (dashed line) and with BQV [|T8 



shadowing and parton energy loss (dashed-dotted line). Parton energy loss is taken to be 



energy dependent, e{ = a s y fi 2 \ a Ei, with /i 2 = lGeV 2 and A a = lfm. 

In Figure (0) we show the dependence of the photon cross section on the form of energy 
loss used, by considering different values of the constant energy loss. For comparison, we 
also show the energy- dependent energy loss case with /i 2 = lGeV 2 and X g = \ q = lfm. It 
is clear that the magnitude of the nuclear cross section is sensitive to the magnitude of the 
energy loss. 

To illustrate the importance of doing the full next-to-leading order calculation of prompt 
photon production at the LHC, we show the nuclear i^-factors, defined as the ratio of 
NLO/LO and as the ratio of NLO/(LO + Br ems.) cross sections in Figure (H). It should 
be clear that the higher order corrections (beyond NLO) are important at low p t and 
that our analysis is only qualitatively reliable there. The dependence of cross section on 
the choice of renormalization, factorization and fragmentation scales is shown in Figure 
(|). Varying scales from 0.5p t to 2.0p t results in theoretical uncertainty of about 30%. 
This indicates the importance of including higher order terms, as does the large K factor 



values presented in Figure (g) . We refer the reader to |TTJ for a complete analysis of scale 
dependence of prompt photon cross section in pp collisions. 

We find that nuclear effects at LHC are significant and should be easily detectable. 
Nuclear shadowing effects are large and need to be better understood. In Figure fllPf) we 
show the rescaled nuclear cross section using the BQV and EKS98 nuclear shadowing. The 
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Figure 6: Prompt photon cross section in the central rapidity region at \/s = 5.5 TeV 
obtained with EKS98 fill shadowing and without parton energy loss (dotted line), with 



parton energy loss (dashed line), with BQV shadowing [18] and no parton energy loss 
(solid line) and with energy loss (dot-dashed line). Parton energy loss is taken to be 

energy dependent, e J a = ot s \l jj?\ a Ea, with /i 2 = lGeV 2 and A a = lfm. 
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Figure 7: Prompt photon cross section in the central rapidity region at y/s = 5.5 TeV 



obtained with EKS98 shadowing |19| and with different values for parton energy loss 
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Figure 8: a) The nuclear K-f actor defined as the ratio of NLO prompt photon cross section 
in heavy-ion collisions, E-jj§-, to the LO cross section and b) the K-factor defined as the 
ratio of NLO to the LO plus bremsstrahlung cross sections. Parton energy loss is taken to 

be energy dependent, e{ = a s J p?\ a E 3 ai with /i 2 = lGeV 2 and A a = lfm. 
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Figure 9: The uncertainty of the photon cross section due to varying factorization, renor- 
malization and fragmentation scales (Q = M = Qp) between 0.5p 4 and 2p t . Parton energy 
loss is taken to be energy dependent, with jj? = lGeV 2 and A a = lfm. 



14 



difference between the two parameterizations of nuclear shadowing is significant at large pt- 
This raises the possibility that one could study the Q 2 dependence of nuclear shadowing 
by measuring the pt spectrum of prompt photons since this difference is independent of the 
form of the energy loss per scattering used. 
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Figure 10: The rescaled prompt photon cross section in the central rapidity region at 
y/s = 5.5 TeV using BQV |H| and EKS98 |T]| shadowing and energy- dependent energy 



loss with fi = lGeV and A a = lfm. 



Current parameterizations of nuclear shadowing in the kinematic region appropriate to 
LHC energies are extrapolations from low energy fixed target data and subject to large 
uncertainties. A precise quantitative knowledge of the nuclear structure functions at the 
small x, large Q 2 kinematic region is crucial. A lepton-nucleus collider such as eRHIC is 
urgently needed. 

Energy loss effects are also large. In Figure (|TTD we show the rescaled nuclear cross 
section for the case of constant energy loss. Whether one can distinguish experimentally 
between different energy loss scenarios will depend on our precise knowledge of nuclear 
structure functions and also on the precision of the prompt photon measurements at RHIC 
and LHC. Distinguishing prompt photons from those coming from decays of pions and 
eta's is notoriously difficult. Also, one expects that the calculation of the ratio of prompt 
photons to pions would reduce some of theoretical uncertainties such as scale dependence 
and intrinsic k t effects, thus making the NLO calculation even more reliable [p3|1 . We intend 
to calculate this ratio in the near future [24]. 
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Figure 11: The rescaled prompt photon cross section in the central rapidity region at 
y/s = 5.5 TeV obtained using constant energy loss and EKS98 |T9] nuclear shadowing 
form. 
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